Delay-induced multiple stochastic resonances on scale-free neuronal networks 
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We study the effects of periodic subtliresliold pacemaker activity and time-delayed coupling on stochastic 
resonance over scale-free neuronal networks. As the two extreme options, we introduce the pacemaker re- 
spectively to the neuron with the highest degree and to one of the neurons with the lowest degree within the 
network, but we also consider the case when all neurons are exposed to the periodic forcing. In the absence 
of delay, we show that an intermediate intensity of noise is able to optimally assist the pacemaker in imposing 
its rhythm on the whole ensemble, irrespective to its placing, thus providing evidences for stochastic resonance 
on the scale-free neuronal networks. Interestingly thereby, if the forcing in form of a periodic pulse train is 
introduced to all neurons forming the network, the stochastic resonance decreases as compared to the case when 
only a single neuron is paced. Moreover, we show that finite delays in coupling can significantly affect the 
stochastic resonance on scale-free neuronal networks. In particular, appropriately tuned delays can induce mul- 
tiple stochastic resonances independently of the placing of the pacemaker, but they can also altogether destroy 
stochastic resonance. Delay-induced multiple stochastic resonances manifest as well-expressed maxima of the 
correlation measure, appearing at every multiple of the pacemaker period. We argue that fine-tuned delays and 
locally active pacemakers are vital for assuring optimal conditions for stochastic resonance on complex neuronal 
networks. 

PACS numbers: 05.45.-a, 05.40.-a, 89.75.Kd 



It is well known that noise can play a constructive role 
in different types of nonlinear dynamical systems, and 
stochastic resonance is perhaps the most prominent exam- 
ple of this fact. The objective of this article is to extend 
the scope of stochastic resonance to complex networks, 
whereby the deterministic periodic input is not only lim- 
ited in its strength but also its outreach. More precisely, 
scale-free neuronal networks are studied on which the sub- 
threshold periodic forcing is introduced only to a single 
neuron of the network, thus acting as a pacemaker. We 
want to determine to what extent the complex scale-free 
topology can aid the pacemaker to entrain the complete 
neuronal ensemble with the help of fine-tuned additive 
noise. Moreover, the new findings are compared with re- 
sults obtained via the more traditional setup where every 
neuron of the network is subjected to a weak periodic forc- 
ing. It is found that scale-free topologies are very efficient 
in propagating noise-supported locaUzed weak rhythmic 
activities. Also, it is found that these paced networks 
are superior to globally forced networks in that stochas- 
tic resonance is better expressed on the former. Impor- 
tantly, since time delays are inherent to the nervous sys- 
tem we take this expUcitly into account via time-delayed 
coupling. We report on the occurrence of delay-induced 
multiple stochastic resonances on scale-free neuronal net- 
works, which appear due to the locking between the delay 
length and the oscillation period of the pacemaker. 



I. INTRODUCTION 

The constructive role of noise in nonlinear dynamics has 
been the subject of intensive research in the past, and it re- 
mains a vibrant topic today. Especially the phenomenon of 
stochastic resonanc e ''^'^ has been studied extensively due to 
its applications in many different fields, ranging from physical 
to social systems. In general, stochastic resonance is charac- 
terized by the optimization of the output signal-to-noise ratio 
in a nonlinear dynamical system following the addition of a 
weak external signal. "^•^•^'^''^ Coherence resonance is a closely 
related phenomenon, which however, refers to the resonant 
response of a dynamical system to pure noise, hence it has 
also been termed as autonomous stochastic resonance i^iiSiii. 
Notably, both stochastic and coherence resonances have been 
reported existing in a wide variety of neuronal models /-^-'°-'^ 
Of particular interest for the present work, due to the usage 
of the same mathematical model of neuronal dynamics, mul- 
tiple resonances in the Rulkov neuronal model'^ have also 
been observed and reported.'^ Effects of correlated noise on 
the dynamics of coupled neurons have also been investigated 
thoroughly, 'iii^ii^iiiiS whereby the importance of different 
temporal and spatial coiTelation lengths has been firmly es- 
tablished. 

Recently, stochastic and coherence resonances on complex 
neuronal networks have attracted increasing attention, and in- 
deed some new features due to the underlying complex in- 
teraction topologies have been revealed. The stochastic res- 
onance on excitable Watts-Strogatz small- world networks - 
via a pacemaker was studied by means of the discrete Rulkov 
mapr^ where it was shown that only for intermediate cou- 
pling strengths is the small-world property able to enhance 
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the stochastic resonance. Moreover, the stochastic reso- 
nance on Newman- Watts networks^' of Hodgkin-Huxley neu- 
rons could be ampHfied via fine-tuning of the small-world 
network structure, depending significantly on the coupling 
strengths among neurons and the driving frequency of the 
pacemaker.^^ Coherence resonance on Watts-Strogatz small- 
world Hodgkin-Huxley neuronal networks has also been in- 
vestigated, and it was found that increasing the randomness 
of the network topology leads to enhancement of temporal 
coherence. Furthermore, Kwon et al?'^ showed that the co- 
herence resonance can be considerably improved by adding 
a small fraction of long-range connections for an intermedi- 
ate coupling strength in a Watts-Strogatz small-world neu- 
ronal network with spatially correlated noise input. Preceding 
these studies were reports on the array-enhanced coherence 
resonance by using an array of coupled FitzHugh-Nagumo 
neurons,ii^i2^ as well as stochastic resonance on small-world 
networks of overdamped bistable oscillators.^* Noise-induced 
phenomena in two-dimensional spatially extended neuronal 
networks have also attracted considerable attention in the past, 
and several studies were devoted to the exploration of possi- 
ble effects of noise. There are some comprehensive reviews 
dedicated to this subject. ^^'-^^ In particular, while spatial co- 
herence resonance was first introduced near pattern forming 
instabilities,-^ it subsequently was reported also in excitable 
neuronal media. Noteworthy, spatiotemporal coherence of 
noise-induced patterns has also been investigated on a regular 
Hodgkin-Huxley neuronal network^^ and it was found that 
the order of the firing rate function could be enhanced as the 
connections amongst neurons became stronger 

In this study, the objective is to extend the scope of stochas- 
tic resonance on complex neuronal networks, particularly 
scale-free networks,-^^ pacemakers, and time-delayed cou- 
pling, thus bringing the setup closer to actual conditions. De- 
lays are inherent to the nervous system because of the fi- 
nite speed at which action potentials propagate across neu- 
ron axons, and due to time lapses occurring in both den- 
dritic and synaptic processing. Notably, it has been sug- 
gested that time delays can facilitate neural synchronization 
and lead to many interesting and even unexpected dynami- 
cal phenomena. ^4.35,36 Moreover, since a power-law distribu- 
tion of the degrees of neurons has been found applicable for 
the coherence among activated voxels using functional mag- 
netic resonance imaging,-'^ and the robustness against simu- 
lated lesions of anatomic cortical networks was also found to 
be very similar to that of a scale-free network,^** our study ad- 
dresses a relevant system setup. We report on the pacemaker- 
driven stochastic resonance in scale-free neuronal networks, 
and compare it with the more classical setup entailing the sub- 
threshold periodic forcing of all neurons constituting the net- 
work. Interestingly, we find that stochastic resonance with 
the locally acting pacemaker is better expressed than the one 
with the globally forced neuronal network. Primarily though, 
we present some non-trivial effects induced by finite delays 
in coupling. In particular, we show that multiple stochastic 
resonances can occur on scale-free neuronal networks if the 
duration of the delay is appropriately tuned. This is primar- 
ily attributed to the emergence of locking between the delay 



length and the oscillation period of the pacemaker. 

The remainder of this paper is organized as follows. In the 
next section the Rulkov map model is reviewed,'^ which will 
be employed to obtain an efficient setup for simulating neu- 
ronal dynamics on scale-free networks. In Section II, the 
time-delayed coupling scheme and the measure for stochastic 
resonance are introduced as well. Main results are presented 
in Section III, whereas the last Section summarizes the new 
findings and discusses their implications. 

n. MATHEMATICAL MODEL AND SETUP 

To effectively simulate the neuronal dynamics on scale- 
free networks, the Rulkov map is utilized,'^ which succinctly 
captures some of the main dynamical features of the com- 
plex continuous-time models yet ensures exceptionally effi- 
cient numerical analysis and processing. The spatiotemporal 
evolution of the studied neuronal network is governed by the 
following iteration equations: 

-^D^e''^ [x^{n-T)-x\n)\, (1) 
j 

= ?/^)(n)-/?a;(^)(n)-7, i = l,...,A/' 

Here, x*^*^ (n) is the membrane potential of the i-th neuron and 
is the variation of ion concentration, representing the 
fast and the slow variables, respectively. The slow temporal 
evolution of y^^^{n) is due to the small values of the positive 
parameters /3 and 7, which within this study are chosen as 
/3 = 7 = 0.001 unless otherwise stated. Moreover, n is the 
discrete-time index, while a is the main parameter determin- 
ing the dynamics of individual neurons constituting the scale- 
free network. If a < 2.0 then all neurons occupy excitable 
fixed (equilibrium) points [x* ~ —^-.V* = — 1 — (a/2)], 
whereas if a > 2.0 then regular bursting oscillations (limit cy- 
cles in the phase space) emerge via a Hopf bifurcationji^ For 
larger a still, these regular oscillations become chaotic3i£ 
Here, we set a = 1.95 and initiate each neuron from fixed 
point initial conditions, so that the additive spatiotemporal 
Gaussian noise ^(''(n), having mean < ^'■'•'(n) >i^n= and 
autocorrelation < S,''^^ in)^^J\h) >— Sij6{n — h), acts as the 
only source of large-amplitude excitations. Furthermore, in 
Eq. (1) f{x) = ^^^2 is a nonlinear function describing the 
essential ingredients of neuronal dynamics, D is the coupling 
strength, the parameter a determines the noise intensity, and 
T is the delay length. The latter three parameters will be in the 
focus of attention within this work, whereas f3 and 7 will be 
varied only occasionally. 

As the base for interactions between neurons the scale- 
free network is employed which is generated via growth and 
preferential attachment as proposed by Barabasi and Albert,^- 
consisting of = 200 vertices in this study. Each vertex cor- 
responds to one neuron, whose dynamics is governed by the 
noise-driven Rulkov map. In Eq. (1), e*'^ = 1 if neuron i is 
coupled to neuron j, and = otherwise. The preferential 
attachment is introduced via the probability 11, which states— 
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that a new vertex will be connected to vertex i depending on 
its degree fcj according to n(fcj) = ki/^- kj. This growth 
and preferential attachment scheme yields a network with an 
average degree kav = ' , and a power-law degree distri- 
bution with the slope of the hne equahng k, —3 on a double- 
logarithmic graph in the present study. Notably, the analytical 
slope of the line is —3 exactly. Such networks having kav = 6 
will be used throughout this work. 

It remains of interest to mathematically introduce the sub- 
threshold periodic pacemaker, which takes the form of a spike 
train defined by 

^(r)(„x Iff' if{nmodt)>{t-w), 
lo, else. 

In Eq. (2) the parameters t, w, g and n are defined as follows: 
t determines the oscillation period of the spike train (pace- 
maker), w is the width of each spike, g defines the amplitude 
of the spikes (the baseline is 0), and n is the discrete time in- 
dex introduced above in Eq. (1). Moreover, the subscript r 
denotes the chosen neuron among all TV = 200 neurons con- 
stituting the excitable (each neuron occupies a fixed point and 
can thus be excited by weak perturbations) scale-free network, 
to which the pacemaker is introduced as an additive term to the 
variable a;*^'"' (n). In our numerical simulations, we choose the 
parameter values t = 700 (unless otherwise stated), w = 50 
and g = 0.015, which ensure that without introducing any 
noise (a = 0) the pacemaker is subthreshold, meaning that by 
itself it cannot induce large-amplitude excitations from any of 
the neurons constituting the network. 

To quantitatively characterize the collective response of the 
neuronal network, we introduce the average membrane poten- 
tial X{n) = X]i=i ('^) the main output to be exam- 
ined further. The correlation of the average membrane poten- 
tial X{n) with the frequency of the pacemaker oj = 27r/t is 
computed via the Fourier coefficients according to 

2 ^* 

Qsin = if^^^i^) sin(wn) (3) 

n=l 



<3cos = ^ 2^ ^l") cos(a;n) (4) 

n=l 

where T = 300 is the number of periods of the pacemaker 
used. Note that the sine and cosine functions in Eqs. (3) 
and (4) have the same frequency uj =- 27r/t as is used for 
pacing the neuronal network [note that t in Eq. (2) is the os- 
cillation period of the spike train], and that thus these equa- 
tions determine the correlation between the output of the net- 
work X(n) and the frequency of the pacemaker. We therefore 
use the Fourier coefficient Q = a/ Q^j^ + Qcos ^s a numer- 
ically effective measure for stochastic resonance, capturing 
succinctly the collective spatiotemporal behavior of the neu- 
ronal network and its correlation with the pacemaker rhythm. 
In general, Q can exhibit a bell-shaped dependence as a key 
parameter (for example a) is varied, indicating the occurrence 




* «e«: vm r$m irtw; 

n 



FIG. 1: Space-time plots obtained at t = and D = 0.006 for 
different a values, equaling: (a) 0.005, (b) 0.008, (c) 0.025, (d) 0.05 
and (e) 0.08. The color profile is linear, black depicting a;^*'(n) = 
—1.7 and white x'-'^ {n) = 0.1, with eight shades of gray in-between. 



of stochastic resonance. Importantly, since the generation of 
scale-free networks has inherent random ingredients, which 
can be additionally amplified by individual vertex (neuron) 
pacing, final results shown below were averaged over 30 in- 
dependent runs for each set of parameter values (wherever ap- 
plicable) to guarantee an appropriate accuracy. 



III. RESULTS 

We start by setting t = and introducing the pacemaker to 
the neuron with the lowest degree kmin within the network, 
thus r = i{kmin)- Space-time plots obtained by varying a are 
presented in Fig. 1, where it can be observed that the excita- 
tory fronts follow the pacemaker rhythm (oscillation period is 
t = 700) only with an intermediate value of a [see Fig. 1(c)]. 
On the other hand, smaller a either largely fail to evoke ex- 
citations [see Fig. l(a)J or the excitatory fronts have defects 
and are not frequent enough [see Fig. l(b)|. Noise intensi- 
ties exceeding the optimal value, however, have an ability of 
initiating excitations on their own (even when the pacemaker 
is not firing), thus again failing to conform to the weakly im- 
posed rhythm [see Fig. 1(d)], or further still, completely over- 
ruling the deterministic dynamics [see Fig. 1(e)]. Presented 
results therefore indicate a classical stochastic resonance sce- 
nario, where an intermediate noise intensity ensures the best 
response of the system to weak external deterministic forc- 
ing. Notably, qualitatively identical results are obtained if the 
pacemaker is introduced to the neuron with the largest degree 
r = i{kmax) within the scale-free network. It is of interest to 
assess these observations quantitatively. 

To describe the pacemaker-driven stochastic resonance 
more precisely, we consider the dependence of Q on cr in Fig. 
2(a) for different options with respect to the placing of the 
pacemaker. It can be observed that, irrespective of whether 
the pacemaker is introduced to the neuron with the minimal 
r = i{kmin) or the maximal degree r = i{kmax), or to all 
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FIG. 2: (a) Dependence of Q on ct at r = and D — 0.006 for different placings of the pacemaker within the scale-free network (see also the 
main text for further details), (b) Dependence of Q on cr and D at r = when the pacemaker is introduced to one of the neurons having the 
lowest degree [r = i(fc„i„)]. 



neurons of the network ["all" in Fig. 2(a)], there exists an 
intermediate optimal noise intensity a for which Q is max- 
imal, thus exhibiting a bell-shaped dependence characteristic 
of the stochastic resonance. Interestingly, one can observe that 
as the pacemaker is introduced to all neurons constituting the 
network, the stochastic resonance decreases rather dramati- 
cally. On the other hand, whether the pacemaker is introduced 
to the neuron with the highest [r = i{kmax)] or to one of 
the neurons with the lowest degree [r = i{kmin)] does not 
notably affect the outlay and magnitude of Q in dependence 
on (T. The placing of an individual pacemaker is mostly ir- 
relevant because in a scale-free network chances are high that 
the low-degree neurons will be connected to one of the main 
high-degree neurons (hubs) with no more than a single direct 
link. Hence now, whether we introduce the pacemaker di- 
rectly to the main hub or one link further away (to any of the 
low-degree neurons) is not vital since the pacemaker rhythm 
does not deteriorate much within nearest or even next-nearest 
neurons. Somewhat more surprising and counterintuitive is 
the fact that the weak periodic forcing of all neurons is less 
efficient in imposing a certain rhythm of excitatory fronts than 
an individual pacemaker We argue that this effect is a conse- 
quence of the essentially non-biased state of neurons consti- 
tuting the scale-free network, where only an individual neuron 
is forced. Such non-biased (not inclined towards an excitation 
by a pacemaker) noisy neurons can respond even to the weak- 
est input from their neighbors, thus synchronizing optimally 
into practically perfect excitatory fronts. On the other hand, 
if all neurons are weakly forced, chances for phase slips are 
much higher since every neuron acts on its own (due to its in- 
dividual forcing), indeed trying to enforce its rhythm to the 
neighbors. Thus, competition between excitations emerges, 
which acts detrimental on the overall synchrony of the net- 
work, in turn decreasing the correlation with the forcing fre- 
quency and leading to the decline of stochastic resonance. In 
fact, the phenomenon shown in Fig. 2(a) is one of the few 
examples where less is actually more in terms of efficiency of 
noise-induced synchrony and coiTelation with weak external 



forcing in neuronal dynamics. 

Before turning to the impact of finite, i.e. non-zero delays 
T, we investigate in Fig. 2(b) stochastic resonance in the tran- 
sition to the strong coupling region (high vales of D). Re- 
sults are presented for r = i{k„iin), but they are qualitatively 
identical also for r ~ i{kmax), and if all the neurons are ex- 
posed to the weak periodic forcing (not shown). Evidently, 
the stochastic resonance phenomenon prevails irrespective of 
D, except that the optimal value of a shifts to slightly higher 
values upon its increase. Thus, it can be concluded that the 
pacemaker-driven stochastic resonance on scale-free neuronal 
networks is a robust phenomenon, occurring largely inde- 
pendently of the particular placing of the pacemaker or the 
strength of the coupling. Moreover, an individually paced 
neuron is actually more effective in warranting ordered ex- 
citatory fronts in accordance with the weakly imposed rhythm 
than the global forcing of the whole network. In what follows, 
therefore, focus is on individually paced noise-driven scale- 
free neuronal networks, where one of the neurons having the 
lowest degree will be chosen as the input for the deterministic 
forcing. 

As when examining the dependence of the neuronal dy- 
namics on (7 in Fig. 1, we start by presenting in Fig. 3 the 
space-time plots obtained for different values of r while keep- 
ing the coupling strength D — 0.006 and the noise intensity 
(T = 0.025 fixed. Results shown in the five panels of Fig. 3 il- 
lustrate the spatiotemporal dynamics of neurons on the studied 
scale-free neuronal network with r ~ i{kmin)- From careful 
visual inspection, an intermittent pattern of regularity and dis- 
order can be inferred upon increasing r. In particular, while 
for T = 700 [panel (b)] and r = 1400 [panel (d)] the excita- 
tory fronts are synchronous and largely obeying to the pace- 
maker rhythm, for r ~ 300 [panel (a)], r — 1000 [panel (c)] 
and T = 1600 [panel (e)] the regularity is either completely 
lost or at least the excitatory fronts become ragged and grad- 
ually loosing synchrony with the imposed frequency. Indeed, 
the delay-induced transitions to ordered spatiotemporal dy- 
namics on scale-free neuronal networks appear intermittently 
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FIG. 3: Space-time plots obtained at D = 0.006 and cr = 0.025 for 
different r values, equalling: (a) 300, (b) 700, (c) 1000, (d) 1400 
and (e) 1600. In all panels the color profile is identical as in Fig. 1. 



FIG. 4: Dependence of cr on r for different combinations of /3 and 7 
(see also main text for further details) when the pacemaker is intro- 
duced to one of the neurons with the lowest deg 
Where applicable, other parameter values are the same as in Fig. 3. 



at roughly integer multiples of the period of the pacemaker, 
equalling t = 700, which corresponds rather accurately to 
the so-called global-resonant oscillation period^ of each indi- 
vidual RuUcov neuron for the currently used parameter values 
(3 — ^ — 0.001 and a = 1.95, which can be extracted from 
the map by calculating the Fourier transform of noise-driven 
oscillations.'*' By setting t largely different from the global- 
resonant oscillation period of individual neurons, however, the 
intermittent outlay presented in Fig. 3 can no longer be ob- 
served. Visual inspection of Fig. 3 thus reveals that regular 
and irregular front propagations appear intermittently as the 
delay is increased, indicating that finite (non-zero) delays in 
coupling among neurons might play a pivotal role in generat- 
ing spatiotemporal order of neuronal activity on scale-free net- 
works in accordance with a weak localized deterministic in- 
put, provided that the latter is adjusted to approximately agree 
with the global-resonant oscillation frequency of the neurons. 

To account for the above visual interpretation quantita- 
tively, we adjust the local dynamics of each neuron (thus far 
we have not varied this) by varying /3 and 7, in turn affect- 
ing the speed of the temporal evolution of yi{n) and con- 
sequently the global-resonant oscillation frequency. At the 
same time, we adjust the oscillation period t of the pace- 
maker correspondingly. In particular, we consider three dif- 
ferent cases; namely f3 = •y ^ 0.0006, /3 — j — 0.001 and 
(3 = J — 0.0016, and change the oscillation period of the 
pacemaker to t — 1200, t — 700 and t — 500, respectively. 
These t values are in good agreement with the global-resonant 
frequency of an individual neuron with the corresponding val- 
ues of (3 and 7. Results presented in Fig. 4 show that, in accor- 
dance with the visual inspection of Fig. 3, multiple resonance 
in Q depending upon the increase of r are obtained by given a 
and D. Following the common terminology, these are termed 
as delay-induced multiple stochastic resonance on scale-free 
neuronal networks. Moreover, it is clear that the particular lo- 
cations of the maxima of Q shift to different values of r as 
/3, 7 and t are varied. Crucially however, it is always so that 
the locking between r and integer multiples of t is preserved. 



Thus, resonances depending on r appear at integer multiples 
of t only if the latter is close to the global-resonant oscillation 
period^ of the individual neurons. On the other hand, val- 
ues of T outside the regions of multiple integers of t impair 
the stochastic resonance significantly, as can be inferred from 
the rather sharp descents of Q towards smaller values, as soon 
as the optimal r's are replaced by other values. We therefore 
conclude that the delay-induced stochastic resonances of neu- 
ronal activity are due to the locking between the delay length 
r and the global-resonant oscillation period of individual neu- 
rons if the latter is close to the oscillation period of the pace- 
maker. This is valid independently of the particular placing of 
the pacemaker, and also for globally paced scale-free neuronal 
networks with different coupling strengths. 

Finally, we present some results where the above-outlined 
dependencies can be observed at a glance. Figures 5(a) and 
(b) show Q in dependence on cr and r for two combinations 
of (3 and 7. In panel (a), where /3 = 7 = 0.001 and t = 700, 
multiple stochastic resonances are clearly eligible as narrow 
white-shaded regions, appearing roughly at integer multiples 
of r = 700 across suitable spans of cr. Practically identical, 
for f3 — "f — 0.0016 and t — 500, stochastic resonances ap- 
pear at integer multiples of r = 500, again strengthening our 
argumentation with respect to the connection with the global- 
resonant oscillation period of individual neurons and the cor- 
responding similarly-set forcing period of the pacemaker It is 
worth emphasizing that regions of optimal r are very narrow, 
especially if compared to the rather broad regions of noise in- 
tensities cr that still warrant reasonably high Q, suggesting that 
fine-tuned delays might be crucial for efficient recognition of 
weak localized external signals. Also notable are some fuzzily 
expressed regions of stochastic resonance at r « 70 in both 
panels of Fig. 5. However, these are most likely due to the fact 
that shorter, if compared to the optimal, delay lengths do not 
influence the neuronal dynamics strong enough to fully pro- 
hibit noise-induced correlations between the pacemaker and 
the neuronal dynamics. This occasionally yields results sim- 
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FIG. 5: Contour plots of Q in dependence on r and a when the pacemaker is introduced to one of the neurons with the lowest degree 
[r = i(fc™i„)] and D = 0.006. Other parameter values are: (a) /? = 7 = 0.001, t = 700 and (b) /? = 7 = 0.0016, t = 500. 



ilar to the case of r = though only at substantially higher 
noise intensities, which are needed to compensate for the dis- 
turbing impact of non-optimal delays. Moreover, we note that 
for other forcing frequencies of the pacemaker we have per- 
formed similar investigation, yet only when the pacemaker 
frequency is close to the global -resonant frequency"" of in- 
dividual neurons forming the scale-free network can multiple 
stochastic resonances be observed at integer multiples of the 
forcing period. 

IV. SUMMARY AND DISCUSSION 

In sum, we have studied stochastic resonance phenomena 
on scale-free neuronal networks in dependence on the noise 
intensity and time-delayed coupling when the pacemaker is 
introduced to different individual neurons or acts as global 
forcing. We found that the stochastic resonance occurs ir- 
respective of the location of the pacemaker and of local or 
global forcing considerations. Yet remarkably, locally intro- 
duced weak pacemakers guarantee better expressed stochastic 
resonance than global forcing. Indeed, it is generally believed 
that global input is not common in real neuronal systems, and 
in fact local inputs are far more likely. In particular, given a 
huge number of neurons, it is unnecessary and even impossi- 
ble to introduce external signals or stimuli to all. Only weak, 
partial and local inputs are reasonable, guaranteing low energy 
consumption and efficiency in large neuronal networks. Fur- 
thermore, by introducing delays to the coupling scheme, we 
observed multiple stochastic resonances upon fine-tuning of 
the delay length, which appear at every multiple of the forcing 
frequency if the latter is close to the global-resonant oscilla- 
tion frequency of individual neurons [note that the determinis- 
tic dynamics is of fixed (equilibrium) point type] 4i More pre- 



cisely, the multiple stochastic resonances appear in an inter- 
mittent fashion as the delay increases, where the intermittency 
is a direct consequence of the on/off locking between the forc- 
ing frequency and the delay length. Thus, we have shown that 
noise and time-delayed coupling play complementary roles in 
ensuring optimal detection of weak localized stimuli in scale- 
free neuronal networks via stochastic resonance. We there- 
fore believe that fine-tuned delays can effectively supplement 
recently identified mechanisms for the enhancement of neu- 
ronal synchronization,"*^'"*^ as well as synchronization"*"* - "*^'"*^'"*^ 
and detection of weak signals^ on complex and scale-free 
networks in general, thereby constituting an important factor 
of interneuronal communication. This argumentation seems 
to be supported also by actual biological data, demonstrat- 
ing that conduction velocities along axons connecting neurons 
vary from 20 to 60 m/s."*^ Real-life delays are thus within the 
range of milliseconds, suggesting that substantially lower or 
higher values may be preclusive for optimal functioning of 
neuronal tissues. We hope our study will be a useful supple- 
ment to the existing body of literature for the function-follow 
form concept,^"'^' as well as the role of complex neuronal net- 
works in generalr^i^ and also serve as a viable source of in- 
formation when striving towards further advances in the field. 
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